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34 Abstract 

35 Both type III effector proteins and non-ribosomal peptide toxins play important 

36 roles for Pseudomonas syringae pathogenicity in host plants, but whether and how 

37 these virulence pathways interact to promote infection remains unclear. Genomic 

38 evidence from one clade of P. syringae suggests a tradeoff between the total number of 

39 type III effector proteins and presence of syringomycin, syringopeptin, and syringolin A 

40 toxins. Here we report the complete genome sequence from P. syringae CC1 557, 

41 which contains the lowest number of known type III effectors to date and has also 

42 acquired genes similar to sequences encoding syringomycin pathways from other 

43 strains. We demonstrate that this strain is pathogenic on Nicotiana benthamiana and 

44 that both the type III secretion system and a new type III effector family, hopBJI, 

45 contribute to virulence. We further demonstrate that virulence activity of HopBJI is 

46 dependent on similar catalytic sites as the E. coli CNF1 toxin. Taken together, our 

47 results provide additional support for a negative correlation between type III effector 

48 repertoires and the potential to produce syringomycin-like toxins while also highlighting 

49 how genomic synteny and bioinformatics can be used to identify and characterize novel 

50 virulence proteins. 
51 

52 Introduction 

53 

54 Pseudomonas syringae, a bacterial phytopathogen of many crop species, utilizes 

55 a diverse arsenal of virulence factors to infect host plants (Hirano and Upper, 2000; 
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56 O'Brien et al., 201 1). Much of the previous research on pathogenicity in P. syringae has 

57 focused on identifying and characterizing virulence factors within strains and on 

58 categorizing their presence across the species (Baltrus et al., 201 1 ; Hwang et al., 2005; 

59 Lindeberg et al., 2012; O'Brien et al., 201 1). Although this accumulation of knowledge 

60 has enabled tests of protein function during infection, many questions remain as to how 

61 virulence pathways interact and evolve at a systems level (Baltrus et al., 201 1 ; Cunnac 

62 et al., 201 1 ; Hwang et al., 2005). In the absence of direct tests, identification of similar 

63 genomic trends across independent lineages provides strong support for the presence 

64 of such interactions over evolutionary time scales. 
65 

66 One of the main contributors to virulence within P. syringae is a type III secretion 

67 system (T3SS) (O'Brien et al., 201 1). The T3SS encodes a structure that translocates 

68 bacterial effector proteins (T3E) into host cells to disrupt host physiological pathways 

69 and enable successful infection (Lindeberg et al., 2012). P. syringae genomes typically 

70 contain about 1 0-40 total T3Es, the exact actions of which depend on the totality of 

71 effector repertoires as well as host genotype (Baltrus et al., 201 1 ; Lindeberg et al., 

72 2012; O'Brien et al., 201 1). Another main contributor to pathogenesis within many P. 

73 syringae strains is non-ribosomal peptide (NRP) toxin pathways. NRP pathways encode 

74 proteins that act as an assembly line, elongating and decorating peptide chains involved 

75 in pathogenicity (Bender et al., 1999). NRP pathways are often not conserved between 

76 closely related strains, but even if present, can differ in their outputs based on 

77 nucleotide polymorphisms (Baltrus et al., 201 1 ; Bender et al., 1 999; Hwang et al., 2005). 
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78 

79 Multiple reports have suggested functional overlap or phenotypic interactions 

80 between T3E and NRP pathways. For instance, phenotypic virulence functions of the 

81 toxin coronatine can be partially restored by the toxin syringolin A as well as T3Es 

82 HopZ1 and AvrB1 (Cui et al., 2010; Jiang et al., 2013; Melotto et al., 2006; Schellenberg 

83 et al., 2010). Furthermore, within MLSA group II P. syringae strains as defined by 

84 Sarkar and Guttman (2004), there exists a negative correlation between presence of 

85 conserved toxin pathways (syringolin A, syringopeptin, and syringomycin) and size of 

86 the T3E repertoire compared to other phylogenetic clades (Baltrus et al., 201 1 ; Figure 

87 S1). In fact, MLSA group II strains possess the lowest number of T3E among analyzed 

88 genomes across the species. This correlation is strengthened by the observation that 

89 strains from pathovar pisi possess some of the largest T3E repertoires within MLSA 

90 group II and have lost all three toxin pathways (Baltrus et al., 2013). Whether this 

91 negative correlation between NRPs and T3Es reflects an unrecognized difference in 

92 disease ecology is yet to be determined, although strains from MLSA group II are 

93 thought to survive as epiphytes much better than other clades within P. syringae (Feil et 

94 al., 2005). 
95 

96 The hypothesis of a phenotypic tradeoff between T3E and NRP pathways would 

97 be bolstered by discovery of an independent clade of P. syringae that has gained similar 

98 toxins to the MLSA group II strains but which also possesses a reduced effector 

99 repertoire. We have analyzed the genomes of a small clade of P. syringae isolated as 



Baltrus MPMI 



4 



Downloaded from http://biorxiv.org/on September 18, 2014 

100 from environmental sources (Morris et al., 2008; Baltrus et al. 2014), and have found 

101 independent evidence supporting this genomic trend (Figure S1). Specifically, compared 

102 to a closely related outgroup, genomes for two strains appear to have lost T3Es while 

103 gaining a pathway which potentially encodes proteins for the production of a 

104 syringomycin-like toxin. We demonstrate that one of these strains, CC1 557, can infect 

105 Nicotiana benthamiana and cause disease in a T3SS-dependent manner. We further 

106 show that virulence of this strain is significantly increased by the presence of one new 

107 T3E effector family, HopBJI , which shares similar structure and catalytic residues with 

108 the CNF1 from Escherichia coli. Therefore, independent evidence suggests that 

109 acquisition NRP pathways correlates strongly with the loss of T3E families and further 

110 strengthens the idea of an ecological or evolutionary tradeoff between these virulence 

111 factors. 
112 

113 Results 

114 

115 A complete genome sequence for P. syringae CC1557 

116 Using a combination of 100bp lllumina paired end reads and longer PacBio reads, we 

117 assembled a complete genome sequence for strain CC1 557 (Morris et al. 2008). This 

118 genome consists of a 5,728,024 bp chromosome and a 53,629 bp plasmid (Genbank 

119 accession number AVEH02000000). According to PGAAP annotation (Angiuoli et al. 

120 2008), this chromosome contains 5001 ORFs while the single plasmid contains an 

121 additional 67 ORFs. 
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122 

123 P. syringae CC1557 can infect and cause disease in N. benthamiana 

124 P. syringae CC1557 was originally isolated from snow, while the closely related strain 

125 CC1466 was originally isolated from asymptomatic Dodecantheon pulchellum, a 

126 perennial herb (Morris et al., 2008). Using syringe inoculation under standard 

127 conditions, we have demonstrated that P. syringae CC1557 can grow to high population 

128 sizes in the apoplast of N. benthamiana after 3 days of infection (Figure 1). This rate of 

129 growth is similar to the well-established pathogenic strain for N. benthamiana, P. 

130 syringae pv. syringae B728a (Figure S2). Furthermore, these large bacterial population 

131 sizes cause disease symptoms as evidenced by the visibility of tissue collapse (Figure 

132 3). High levels of growth and tissue collapse are both eliminated by deleting the hrcC 

133 gene (Figures 1 and 3), which encodes a main structural protein for the type III secretion 

134 pilus. Therefore, CC1557 virulence in N. benthamiana under laboratory conditions is 

135 T3SS dependent. 
136 

137 The genomes of P. syringae CC1557 and CC1466 encode a reduced type III 

138 effector repertoire. We used bioinformatic methods as per Baltrus et al. 201 1 , to 

139 search the complete and draft genomes of P. syringae strains CC1557, CC1466, and 

140 the closely related strain CC1583, for known type III effector families (Baltrus et al. 

141 2014). We found a total of four T3E shared across these three genomes, with eight T3E 

142 found only within CC1583 (Table 2). Therefore, together with HopBJI (see below), the 

143 genomes of CC1466 and CC1557 appear to encode a total of 5 potential T3E from 
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144 known protein families. This stands in contrast to the immediate outgroup strain 

145 CC1583, which encodes a total of 12 potential T3E from known protein families. 

146 

147 P. syringae CC1466 and CC1557 have horizontally acquired pathways for 

148 production of non-ribosomal peptides that resemble syringopeptin and 

149 syringomycin. We used tBLASTn searches of the complete genome sequence of 

150 CC1 557 in order to identify potential NRP toxin pathways as per Baltrus et al. 201 1 . 

151 This genome contains highly similar full-length tBLASTn matches to a variety of proteins 

152 involved in syringomycin synthesis: SyrB (93% Protein Sequence Identity), SyrC (99%), 

153 and SyrP (93%), data not shown. Further investigation of all complete genomes for P. 

154 syringae demonstrates that the region containing putative NRP-related loci in CC1 557 is 

155 found in approximately the same genomic context on the chromosome as the 

156 syringomycin-syringopeptin pathways in P. syringae pv. syringae B728a (Figure S3). 

157 Moreover, the genomic context surrounding these toxin pathways is conserved 

158 throughout P. syringae (Figure S3). Although regions of high similarity to the 

159 syringopeptin pathway do not appear to be present within CC1557, analysis of gene 

160 annotations suggests that a second NRP pathway occurs immediately downstream of 

161 the syringomycin-like region on the chromosome of CC1 557 (data not shown). While it 

162 is difficult to identify the NRP product of these CC1 557 pathways by sequence alone, 

163 draft genomes sequences suggest that this region is present in CC1466 but not 

164 CC1583. The most parsimonious explanation for this phylogenetic signal is that an 
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165 immediate ancestor of CC1466 and CC1557 acquired this NRP region by horizontal 

166 transfer after divergence from the CC1 583 lineage. 

167 

168 A new type III effector protein family is present in the genomes of both P. 

169 syringae CC1466 and CC1557. Similarity searches of the P. syringae genomes 

170 demonstrated that hopM1, which encodes an effector protein largely conserved 

171 throughout P. syringae, had likely been lost within strains CC1466 and CC1557. Further 

172 investigation demonstrated that an unknown protein open reading frame had replaced 

173 hopM1 within these genomes (Figure 2). In order to test if this unknown protein is 

174 translocated into plant cells, we cloned the region encoding the predicted ORF into the 

175 translocation test vector pBAV1 78 and tested for delivery of AvrRpt2 as a fusion protein. 

176 We found that the resulting HopBJ1-AvrRpt2 fusion protein triggered cell death in both 

177 Arabidopsis thaliana Col-0 wild type and rps2 lines, indicating that HopBJI itself was 

178 capable of triggering cell death in the absence of recognition of AvrRpt2 (Table 3, Figure 

179 4B). In order to verify that HopBJI triggered cell death without AvrRpt2 we delivered 

180 HopBJI :HA into a variety of plant genotypes using Pto DC3000 D28E (Cunnac et al., 

181 201 1 ). Delivery of a HopBJI :HA can cause tissue collapse in a variety of Arabidopsis 

182 accessions (Table 3). All genotypes tested displayed strong cell death by either 24 or 

183 48hr post inoculation. 
184 

185 HopBJI induced cell death is not dependent on the known R-gene pathways as 

186 tissue collapse still occurs in edsl, pad4 or rarl mutants (Table 3). HopBJI -induced cell 
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187 death also occurs in plants expressing the salicylic acid degrading enzyme NahG (Table 

188 3). However, this rapid tissue collapse fails to occur out of a hrcC- background of Pto 

189 DC3000, indicating HopBJI secretion is dependent on a functioning T3SS (data not 

190 shown). The genomic region immediately upstream of hopBJI possesses a canonical 

191 hrp-box sequence (Fig. S4), and so regulation of this gene very likely takes place 

192 through the action of the sigma factor HrpL (Mucyn et al., 2014). We therefore conclude 

193 that this novel ORF encodes a new putative T3E family, HopBJ, and that this family 

194 possesses broad cytotoxic capabilities. 
195 

196 Deletion of HopBJI from P. syringae CC1557 causes loss of disease symptoms 

197 and lowers bacterial populations in plants. We have created a deletion of hopBJI in 

198 P. syringae strain CC1557, and have performed growth curves in planta using syringe 

199 infiltration. Deletion of hopBJI causes a repeatable loss of growth in planta after two 

200 days of infection (Figure 3A). We further find that deletion of hopBJI leads to a loss of 

201 the tissue collapse phenotype after 3 days (Figure 3B). We have been able to 

202 complement both of these phenotypes in planta through expression of hopBJI using its 

203 native promoter in trans (Figure 3). 
204 

205 HopBJI resembles cytotoxic genes present within other species. Currently, the 

206 closest BLASTp hits for HopBJI are hypothetical proteins from Serratia marcescens 

207 (34% identity, YP_007346200. 1 ) and Hahella chejuensis (30% identity, YP_4351 89. 1 ) . 

208 However, modeling of protein structure using the Phyre2 web server demonstrates that 
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209 HopBJI displays limited similarity to the E. coli CNF1 toxin (Figure 4A). Surprisingly, 

210 amino acids critical for the deamidase functions of CNF1 appear to be conserved in 

211 HopBJI (Figure 4A) as well as in similar proteins from S. marcescens and H. chejuensis 

212 (Buetow et. al., 2001 ; data not shown). This modeling of protein structure further 

213 suggests that C174 and H192 are a functional catalytic dyad within HopBJI . 
214 

215 Predicted catalytic residues of HopBJI are required for effector function. To test 

216 for functionality of predicted catalytic residues for HopBJI , we measured the effect of 

217 putative catalytic dead mutants (C174S, H192A) in planta. Although mutant versions of 

218 HopBJI are successfully translocated by strain Pto DC3000 D28E into Arabidopsis, only 

219 the wild type version causes tissue collapse in rps2 101c plants (Figure 4B). Moreover, 

220 both catalytic dead mutants are unable to complement virulence defects of a hopBJI 

221 deletion in CC1557 in N. benthamiana (Figures 5C and 5E). Both mutants, unlike the 

222 wild type version, also fail to lower the growth of PtoDC3000 in Arabidopsis (Figure 4D) 

223 and fail to cause tissue collapse when transiently expressed in N. benthamiana. 

224 Therefore, as with similar residues in E. coli CNF1 (Buetow et. al., 2001), we conclude 

225 that both C174 and H192 are required for full enzymatic function of HopBJI . 
226 

227 Discussion 

228 

229 Recent advances in genome sequencing have spurred great interest in 

230 conservation and diversification of virulence pathways across Pseudomonas syringae 



Baltrus MPMI 



10 



Downloaded from http://biorxiv.org/on September 18, 2014 

231 (O'Brien et al., 201 1 ; Studholme, 201 1). While genome gazing often leads to the 

232 identification of trends across a species, single observations must be treated with 

233 skepticism in the absence of additional data. One particularly interesting, yet relatively 

234 little understood, pattern across genomes involves a negative correlation in MLSA group 

235 II between size of T3E repertoire and presence of NRP pathways (Baltrus et al., 201 1 ; 

236 Figure S1 ). Importantly, this trend does not appear to be a result of sampling bias or 

237 failure to discover novel T3E because strains within this clade have been thoroughly 

238 screened at both genetic and transcriptome level (Baltrus et al., 201 1 ; Mucyn et al., 

239 2014). Indeed the entire hrpL regulon, the major transcriptional regulator for all P. 

240 syringae T3E, is reduced within MLSA group II compared to other clades (Mucyn et al., 

241 2014). 
242 

243 Strains CC1 466 and CC1 557 diverge early in the phylogeny of P. syringae sensu 

244 latu, and have not been isolated from diseased plants in nature (Morris et al., 2008). We 

245 have shown that CC1 557 can act as a pathogen of N. benthamiana and, as with all 

246 other phytopathogenic P. syringae strains, that virulence requires a functional T3SS 

247 (Figure 1 ). We have found that both of these genomes contain few T3Es, with only 5 

248 known loci present and shared by both (Table 2). Genomic comparisons between these 

249 two strains have further demonstrated that both contain an NRP toxin pathway very 

250 similar to those encoding syringomycin. As is the nature of NRP toxins even a single 

251 amino acid change can have dramatic consequences on the final toxin product 

252 chemistry (Bender et al., 1 999), but it is important to note that these many of the 
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253 predicted ORFs in these regions are highly similar (>90% protein ID) to syringomycin 

254 ORFs present within MLSA group II strains. Comparison to closely related outgroup 

255 strains demonstrates that acquisition of this NRP pathway is evolutionarily correlated 

256 with loss of T3E families. Indeed, the most closely related sequenced outgroup strain, 

257 CC1583, lacks this NRP pathway and contains twice the number of known T3E loci 

258 (Table 2). As such, we interpret these patterns as independent evidence suggesting an 

259 evolutionary or ecological tradeoff between presence of syringomycin and number of 

260 T3Es within the strain. We also note that there appears to be an additional, 

261 uncharacterized NRP pathway downstream of the syringomycin-like region within 

262 CC1 557 (data not shown). 
263 

264 What selective pressure could underlie such a negative correlation between T3E 

265 and NRP toxins? It is possible that the functions of these NRP toxins are redundant with 

266 a specific suite of T3Es. In this case, T3E loss can take place because there is no 

267 longer positive selection to counter selection against evolutionary (such as recognition 

268 by the plant immune system) or physiological costs. For MLSA clade II and CC1557, 

269 acquisition of syringomycin and potentially additional NRP pathways could render the 

270 functions of some T3E obsolete. If this is the case, we further predict that the potential 

271 second NRP toxin pathway (downstream of the syringomycin-like pathway in CC1557) 

272 could have similar functions as syringopeptin. It is also possible that disease ecologies 

273 of most MLSA group II strains, as well as CC1466 and CC1557, differ compared to 

274 other focal P. syringae strains. We have a rudimentary understanding of differences in 
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275 disease ecology throughout P. syringae pathogens, although it is recognized that 

276 differences do exist across the phylogeny (Clarke et al., 2010). Specifically, one clade of 

277 strains from MLSA group II has exchanged the canonical P. syringae T3SS with a 

278 second divergent T3SS at an independent genomic location while still maintaining toxin 

279 pathways. While these strains can grow in planta, it is unknown whether they can cause 

280 disease (Clarke et al. 2010). Furthermore, infection of woody hosts has arisen multiple 

281 independent times on multiple hosts yet such strains are still dependent on a functioning 

282 T3SS (O'Brien et al., 201 1 ; Ramos et al., 201 2). Strains that contain syringopeptin and 

283 syringomycin potentially cause disease on a wider host range than other groups, which 

284 may be enabled by generality of NRP toxins (Baltrus et al., 201 1 ; Hwang et al., 2005; 

285 Quigley and Gross, 1994). Alternatively, these strains may survive and persist across a 

286 variety of environments and host plants, but may only rarely reach high enough 

287 population densities to cause disease. MLSA group II strains appear to survive better as 

288 epiphytes compared to other groups, and presence of NRP toxins may tie into this 

289 strategy (Feil et al., 2005). Moreover, such differences in disease ecology may be 

290 difficult to precisely measure because they may only be visible under natural conditions 

291 of infection and dispersal. 
292 

293 It is possible that there are a number of undescribed T3E families within CC1466 

294 and CC1 557, and to this point we have successfully used comparisons of genome 

295 synteny to identify HopBJI . hopBJI is present within the conserved effector locus within 

296 these strains, and as such is proximate to both the structural genes for the T3SS and 
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297 the T3E avrE. While it is difficult to precisely map out the evolutionary history of this 

298 region, it does appear that hopBJI has been recombined into this locus in place of 

299 hopM1 (Figure 2). In some respects, recombination of hopM1 from this locus is 

300 analogous to a situation witnessed in MLSA group I P. syringae strains where hopM1 

301 alleles have been cleanly swapped through recombination (Baltrus et al., 201 1). HopM1 

302 is an effector whose presence (although not necessarily in functional form) is conserved 

303 throughout P. syringae sensu strictu, and has been shown to act redundantly with AvrE 

304 during infection of plant hosts (Badel et al., 2006; Kvitko et al., 2009). Moreover, in 

305 contrast to frameshift or nonsense mutations disrupting the coding sequence of hopM1 

306 in P. syringae, presence of this gene is polymorphic within P. viridiflava. While some 

307 strains maintain an intact version, others lack hopM1 completely (Araki et al.,2006; 

308 Bartoli et. al, 2014). It is currently unknown whether HopBJI can act redundantly with 

309 AvrE on specific hosts, although that our single deletion mutant shows a virulence 

310 phenotype on N. benthamiana speaks against this possibility. Moreover, since HopBJI 

311 appears to be responsible for a significant portion the growth of CC1 557 in planta, it will 

312 be interesting to see how this effector behaves within different strain backgrounds on 

313 different hosts. 
314 

315 HopBJI itself displays protein structure similarity to CNF1 toxin found within 

316 pathogenic E. coli strains. CNF1 functions by causing deamidation of a glutamine 

317 residue, crucial for GTP hydrolysis for small GTPases of the Rho family, thereby leading 

318 to constitutive activation and actin disruption (Lemonnier et al., 2007). Furthermore, 
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319 changing either of these amino acids within HopBJI eliminates virulence activity of 

320 CC1 557 in planta, as well as during transient expression by Agrobacterium, and 

321 eliminates avirulence activity when delivered from PtoDC3000 in Arabidopsis (Figure 4). 

322 Potential functional parallels between CNF1 and HopBJI are particularly interesting 

323 because Rho GTPases are known to regulate cytoskeletal dynamics in plants (Mucha et 

324 al., 201 1), actin has a demonstrated role in basal defense (Henty-Ridilla et al., 2013), 

325 and the P. syringae effector HopZ1 a has been shown to target tubulin in order to 

326 promote virulence (Lee et al., 2012). If HopBJI does indeed target Rho family GTPases 

327 to disrupt the cytoskeleton, it would represent a striking example of molecular 

328 convergence across plant and animal pathogens. 
329 

330 The pathology of HopBJI on Arabidopsis, when delivered from strain Pto 

331 DC3000, suggests that HopBJI could be acting as a general toxin because cell death 

332 appears regardless of accession and is not dependent on typical R-gene related host 

333 defense pathways (Table 3). Moreover, HopBJI does act as an avirulence factor within 

334 PtoDC3000 because this cell death limits the growth of this strain in Arabidopsis (Figure 

335 4D). In some ways, this may be similar to the functions of AvrE, which was identified as 

336 causing cell death in soybean (Kobayashi et al., 1989). It is tempting to think that both 

337 AvrE and HopBJI ultimately lead to the same physiological outcomes during infection, 

338 but so far little is known about processes targeted by AvrE. 
339 
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340 What advantages HopBJI provides to CC1557 in nature remains unknown, as it 

341 is currently unclear whether this strain is capable of causing disease in any organism 

342 outside of the lab environment or syringe infiltration. Indeed, both CC1557 (snow) and 

343 CC1466 (asymptomatic plants) were originally isolated in order to characterize 

344 environmental diversity of P. syringae strains outside of crop disease. It is certainly 

345 possible that CC1 557 and related strains persist at low levels across a range of plant 

346 hosts and habitats, using a T3SS and toxins to survive epiphytically but never reaching 

347 large enough population sizes to cause disease. It is also possible that CC1557 is 

348 naturally pathogenic within a suite of under-studied plant hosts or only under specific 

349 environmental conditions. However, given the disease symptoms measured under 

350 laboratory conditions, it is clear that CC1 557 has the potential to cause rapid cell death 

351 in plant leaves relative to other P. syringae strains and that a significant amount of this 

352 activity is dependent on HopBJI . These results suggest an intriguing possibility that the 

353 disease ecology of CC1 557 is fundamentally different (for instance, more necrotrophic) 

354 than other commonly studied P. syringae strains. 
355 

356 Herein we have analyzed the genomes of two P. syringae strains isolated from 

357 environmental sources. These genomes contain a paucity of known T3E, even though 

358 virulence of CC1557 on N. benthamiana is dependent on a functioning T3SS. That 

359 these strains have also acquired a NRP pathway similar to syringomycin represents an 

360 independent evolutionary example that supports a negative correlation between NRP 

361 pathways and T3E repertoires. Furthermore, this data set highlights how genomic scans 
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362 can lead to insights into pathogenesis while also demonstrating the power of genetic 

363 diversity to uncover genome-wide changes in virulence gene architecture. 

364 

365 Methods 

366 

367 Plasmids, bacterial isolates, and growth conditions. 

368 

369 All bacterial strains and plasmids used or created are listed in table 1 . Typically, 

370 P. syringa e and Agrobacteri um isolates were grown at 27°C on KB media using 

371 50ng/mL rifampicin. When necessary, cultures of P. syringa e, Agrobacteri um and E. 

372 coli were supplemented with antibiotics or sugars in the following concentrations: 

373 10ng/ml_ tetracycline, 50ng/ml_ kanamycin, 50ug/ml spectinomycin, 25ng/mL 

374 gentamycin (50ug/mL for Agrobacteri um), and 5% sucrose. 
375 

376 All clones were created by first amplifying target sequences using Pfx 

377 polymerase (Invitrogen), followed by recombination of these fragments into the entry 

378 vector pDONR207 using BP clonase (Invitrogen). Site-directed mutagenesis was 

379 performed using SLIM amplification ( Chiu et al., 2004}. All ORF (without a hrp-box 

380 promoter) and gene (including the promoter) sequences were confirmed by Sanger 

381 sequencing of these pDONR207 clones. Clones in entry vectors were recombined into 

382 destination vectors of interest using LR clonase (Invitrogen). 
383 
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384 Genome Sequences and Searches 

385 

386 Draft genome assemblies for CC1 466, CC1 557,and CC1 583 are publicly 

387 available through Genbank (accession numbers AVEM00000000, AVEHOOOOOOOO, and 

388 AVEGOOOOOOOO respectively; Baltrus et al., 2014). One PacBio SMRT cell yielded 

389 37,509 reads for a total of 268, 1 22,626 nucleotides after filtering for quality. To piece 

390 together the complete genome sequence for CC1 557, PacBio reads were first 

391 assembled using the HGAP software (Chin et al., 2013), yielding two contigs total with 

392 no scaffolding gaps. Contigs from the previous assembly (using only 100bp lllumina 

393 paired end reads, Baltrus et al. 2014) were then overlayed onto this HGAP assembly. 

394 When there was coverage by contigs from the lllumina assembly, the previous 

395 assembly sequence was chosen as the final sequence, but when there was no 

396 coverage the PacBio assembly was chosen. In cases where multiple divergent lllumina 

397 contigs assembled to the same region in the HGAP assembly, the lllumina sequence 

398 which matched the PacBio assembly was chosen. Since the second contig from the 

399 HGAP assembly possesses numerous plasmid related genes, we are confident that this 

400 does actually represent a plasmid present within this strain. Both the chromosome and 

401 plasmid were annotated using PGAAP (Angiuoli et al., 2008). 
402 

403 T3Es and NRP pathways were identified as per Baltrus et al. (201 1). Briefly, draft 

404 genome assemblies were queried using tBLASTn with queries consisting of protein 

405 sequences for known type III effector families or for key loci within toxin pathways. Each 
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406 BLAST hit was validated by hand for copy number, identity, and completeness. All 

407 BLAST results were visually inspected to make sure that each sequence displayed high 

408 similarity to only one region in the assembly, that this similar region was not part of a 

409 larger ORF, and that the length of this region was greater than 40% of the length of the 

410 original query sequence. 
411 

412 The Phyre2 server ( http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index ) 

413 was used to carry out protein structure similarity searches (Kelley and Sternberg, 2009). 

414 In brief, this server performs an automated search of both the protein sequence and 

415 predicted folds for proteins of interest (in this case HopBJI ) across multiple databases. 

416 Within this pipeline, iterative searches such as PSI-Blast enable identification of 

417 distantly similar sequence matches. Once a profile is constructed for the sequence of 

418 interest, it is compared against a database of known protein structures and returns the 

419 best matches as a prediction of protein folds within the query. 
420 

421 Generation of bacterial mutants 
422 

423 Bacterial mutants were generated as per Baltrus et al. 2012. Regions (>500bp) 

424 upstream and downstream of the target genes were PCR amplified separately and then 

425 combined into one fragment by overlap extension PCR. The bridge PCR amplicon was 

426 then cloned into pDONR207, and further moved into pMTN1907 using LR clonase. 

427 Once mated into P. syringae, single recombination of a homologous region upstream or 
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428 downstream of the target region and subsequent selection on sucrose allows for 

429 screening of clean deletions. Mutants were confirmed by phenotyping for sucrose 

430 resistance, tetracycline sensitivity, PCR amplification of the deletion, and failure of PCR 

431 to amplify regions within the deletion. 

432 

433 Negative Correlation Between Syringomycin-like Toxin Pathways and Type III Effectors 

434 

435 We calculated statistical significance for negative correlation between toxin 

436 pathways and effectors using data on the number of full length T3Es per strain from 

437 Baltrus et al. 201 1 as well as data from Table 2 (Figure S2). To minimize bias due to 

438 phylogenetic relationships, we chose > 5 diverse strains from MLSA groups 1 , 2, and 3 

439 as well as CC1557 and CC1583. We then performed a Wilcoxian rank sum test to 

440 compare the number of full-length effectors between genomes which either contain or 

441 lack the genetic capacity for production of a syringomycin-like toxin. 
442 
443 
444 

445 Nicotiana benthamiana plants were grown for 4-6 weeks on a long day cycle (1 6 

446 hours light/8 hours dark). Plants were removed from the growth chamber and 

447 allowed to acclimate to ambient laboratory conditions for a period of 3-5 days prior to 



448 



In planta Growth Curve s, Cell-death assays and Translocation Tests 



infiltration (Arizona) . Bacteria were cultured over night on KB amended with proper 



449 antibiotics. Bacterial cells were washed 1x in and resuspended to an OD 6 oo of 
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450 0.002 in sterile 10 mM MgCI 2 yielding a final inoculum density of approximately 1x10 5 

451 (North Carolina) or 1x10 6 CFU/mL (Arizona) . Cell suspensions were then syringe 

452 infiltrated into the abaxial side of N. benthamiana leaves. Populations were recovered 

453 after 2 (North Carolina) or 3 (Arizona) days of growth using a corkborer unless 

454 otherwise noted. Leaf disks were disrupted using glass beads and a bead-beater device 

455 and populations were enumerated by dilution plating onto KB amended with appropriate 

456 antibiotics. Arabidopsis growth assays were performed similarly, with the exception that 

457 plants were grown in walk-in rooms maintained at on a short day cycle (9 hours light at 

458 21 C and 15 hours dark at 18C). 
459 

460 Cell death induced by hopBJI was assayed after transient expression of 35S- 

461 driven hopBJI :YFP in N. benthamiana. Bacteria were grown to saturation overnight in 

462 2xYT media and then diluted to an OD goo of 0.1 . Cultures were syringe injected into fully 

463 expanded leaves of 5-6 week old plants and cell death was visualized 24hr post- 

464 innoculation. Accumulation of WT and mutant HopBJI was verified with standard 

465 Western blotting of lysates from injected leaf cores using anti-GFP (Roche) 
466 

467 HopBJI was tested for the ability to secrete the active C-terminal fragment of 

468 AvrRpt2 to cause a hypersensitive response (HR) in Arabidopsis accession Col-0 

469 (Guttman et al., 2002). For this test, promoterless hopBJI was cloned into plasmid 

470 pBAV178 (Vinatzer et al., 2006) so that expression in strain Pto DC3000 D28E was 

471 dependent on a tet promoter. This construct was also placed into a hrcC- mutant of Pto 
All DC3000 (Baltrus et al., 2012), in order to test the requirement of a functioning TTSS for 
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473 tissue collapse. All HR tests were performed on ~5 week old plants and utilized a 

474 bacterial density of OD 6 oo 0.05. Tissue collapse was measured at either 24 or 48 hours 

475 post infection and confirmed by at least four independent tests. For translocation of 

476 hopBJI and hopBJI site-directed mutants, similar assays were done with native 

477 promoter clones in pJC532 delivered from the D28E effector-deleted strain of Pto 

478 DC3000. ( Chang et al., 2005) 
479 

480 Phylogenetic Methods 

481 

482 Phylogenetic analyses were performed on concatenated MLSA loci as described 

483 previously using concatenated fragments from 7 MLSA genes across a diverse array of 

484 isolates from P. syringae (Baltrus et al, 201 1 ) as well as AvrE and HrpW protein 

485 sequences. MrBayes was used to perform Bayesian phylogenetic analyses with flat 

486 priors, a burn-in period of 250,000 generations, and convergence after 1 ,000,000 total 

487 generations (Ronquist et al., 2012). 
488 



489 

490 Angiuoi, S.V., Gussman, A., Klimke, W., Cochrane, G., Field, D., Garrity, G., Kodira, 

491 CD., Kyrpides, N., Madupu, R., Markowitz, V., et al. (2008). Toward an online 

492 repository of Standard Operating Procedures (SOPs) for (meta)genomic 

493 Annotation. OMICS, 12, 137-41 

494 Araki, H., Tian, D., Goss, E.M., Jakob, K., Halldorsdottir, S.S., Kreitman, M., and 

495 Bergelson, J. (2006) Presence/absence polymorphism for alternative 

496 pathogenicity islands in Pseudomonas viridiflava, a pathogen of Arabidopsis. 

497 Proc. Natl. Acad. Sci USA 103, 5887-5892. 

498 Badel, J.L., Shimizu, R., Oh, H.-S., and Collmer, A. (2006). A Pseudomonas syringae 

499 pv. tomato avrE1/hopM1 mutant is severely reduced in growth and lesion 

500 formation in tomato. Mol Plant Microbe Interact 19, 99-1 1 1 . 



Baltrus MPMI 



22 



Downloaded from http://biorxiv.org/on September 18, 2014 



501 Baltrus, D.A., Nishimura, M.T., Romanchuk, A., Chang, J.H., Mukhtar, M.S., Cherkis, K., 

502 Roach, J., Grant, S.R., Jones, CD., and Dangl, J.L. (2011). Dynamic Evolution of 

503 Pathogenicity Revealed by Sequencing and Comparative Genomics of 19 

504 Pseudomonas syringae Isolates. PLoS Pathog 7, e1 0021 32. 

505 Baltrus, D. A., Nishimura, M. T., Dougherty, K. M., Biswas, S., Mukhtar, M. S., Vicente, 

506 J., et at. (201 2) The molecular basis of host specialization in bean pathovars of 

507 Pseudomonas syringae. Mol Plant Microbe Interact 25, 877-888. 

508 Baltrus, D.A., Dougherty, K., Beckstrom-Sternberg, S.M., Beckstrom-Sternberg, J.S., 

509 and Foster, J.T. (2013) Incongruence between MLSA and whole genome 

510 phylogenies: Pseudomonas syringae as a cautionary tale. Mol Plant Pathol, doi: 

511 10.1111/mpp.12103 

512 Baltrus, D.A., Yourstone, S., Lind, A., Guilbaud, C, Sands, D.C., Jones, CD., Morris, 

513 C.E., Dangl, J.L. (2014). Draft Genomes for a phylogenetically diverse suite of 

514 Pseudomonas syringae strains from multiple source populations. Genome Announc 

515 23, e01 195-13 doi: 10.1 128/genomeA.01 195-13 

516 Bartoli, C, Berge, O., Monteil, C.L., Guilbaud, C, Balestra, G.M., Varvaro, L, Jones, C, 

517 Dangl, J.L., Sands, D.C, and Morris, CE. (2014) The Pseudomonas viridiflava 

518 phylogroups in the P. syringae species complex are characterized by genetic 

519 variability and phenotypic plasticity of pathogenicity-related traits. Environ Microbiol 

520 doi: 10.1111/1462-2920.12433. 

521 Bender, C.L., Alarcon-Chaidez, F., and Gross, D.C. (1999). Pseudomonas syringae 

522 phytotoxins: mode of action, regulation, and biosynthesis by peptide and polyketide 

523 synthetases. Microbiol Mol Biol Rev 63, 266-292. 

524 Buetow, L, Flatau, G., Chiu, K., Boquet, P., and Ghosh, P. (2001). Structure of the Rho- 

525 activating domain of Escherichia coli cytotoxic necrotizing factor 1 . Nat Struc Biol 8, 

526 584-588 

527 Chang, J.H., Urbach, J.M., Law, T.F., Arnold, L.W., Hu, A., Gombar, S., Grant, S.R., and 

528 Dangl, J.L. (2005). A high-throughput, near saturating screen for type III effector 

529 genes from Pseudomonas syringae. Proc Natl Acad Sci USA 1 02, 2549-2554 

530 Chin, C.-S., Alexander, D.H., Marks, P., Klammer, A.A., Drake, J., Heiner, C, Clum, A., 

531 Copeland, A., Huddleston, J., Eichler, E.E., et al. (2013). Nonhybrid, finished 

532 microbial genome assemblies from long-read SMRT sequencing data. Nat Meth 

533 10, 563-569. 

534 Chiu, J., March, P.E., Lee, R., Tillett, D. (2004) Site-directed, Ligase-independent 

535 mutagenesis (SLIM): a single-tube methodology approaching 100% efficiency in 

536 4h. Nuc Acids Res 32, e1 74 



Baltrus MPMI 



23 



Downloaded from http://biorxiv.org/on September 18, 2014 



537 Clarke, C.R., Cai, R., Studholme, D.J., Guttman, D.S., and Vinatzer, B.A. (2010). 

538 Pseudomonas syringae strains naturally lacking the classical P. syringae hrp/hrc 

539 Locus are common leaf colonizers equipped with an atypical type III secretion 

540 system. Mol Plant Microbe Interact 23, 198-210. 

541 Cui, H., Wang, Y., Xue, L, Chu, J., Yan, C, Fu, J., Chen, M., Innes, R.W., and Zhou, 

542 J.-M. (2010). Pseudomonas syringae Effector Protein AvrB Perturbs Arabidopsis 

543 Hormone Signalingby Activating MAP Kinase 4. Cell Host Microbe 7, 164-175. 

544 Cunnac, S., Chakravarthy, S., Kvitko, B.H., Russell, A.B., Martin, G.B., and Collmer, A. 

545 (201 1). Genetic disassembly and combinatorial reassembly identify a minimal 

546 functional repertoire of type III effectors in Pseudomonas syringae. Proc Natl 

547 Acad Sci USA 108, 2975-2980. 

548 Feil, H., Feil, W.S., Chain, P., Larimer, F., DiBartolo, G., Copeland, A., Lykidis, A., 

549 Trong, S., Nolan, M., Goltsman, E., et al. (2005). Comparison of the 

550 complete genome sequences of Pseudomonas syringae pv. syringae B728a 

551 and pv. tomato DC3000. Proc Natl Acad Sci USA 102, 1 1064-1 1069. 

552 Guttman, D.S., Vinatzer, B.A., Sarkar, S.F., Ranall, M.V., Kettler, G., and Greenberg, 

553 J.T. (2002). A functional screen for the type III (Hrp) secretome of the plant 

554 pathogen Pseudomonas syringae. Science 295, 1 722-1 726. 

555 Henty-Ridilla, J.L, Shimono, M., Li, J., Chang, J.H., Day, B., and Staiger, C.J. (2013) 

556 The plant actin cytoskeleton responds to signals from microbe-associated 

557 molecular patterns. PLoS Pathog 9, e1 003290 

558 Hirano, S., and Upper, C. (2000). Bacteria in the leaf ecosystem with emphasis on 

559 Pseudomonas syringae— a pathogen, ice nucleus, and epiphyte. Microbiol Mol 

560 Biol Rev 64, 624. 

561 Hwang, M.S.H., Morgan, R.L., Sarkar, S.F., Wang, P.W., and Guttman, D.S. (2005). 

562 Phylogenetic characterization of virulence and resistance phenotypes of 

563 Pseudomonas syringae. Appl Environ Microbiol 71, 5182-5191. 

564 Jiang, S., Yao, J., Ma, K.-W., Zhou, H., Song, J., He, S.Y., and Ma, W. (2013). Bacterial 

565 Effector Activates Jasmonate Signaling by Directly Targeting JAZ Transcriptional 

566 Repressors. PLoS Pathog 9, e1 00371 5. 

567 Kelley, L.A., and Sternberg, M.J.E. (2009). Protein structure prediction on the Web: a 

568 case study using the Phyre server. Nat Protoc 4, 363-371 . 

569 Kobayashi, D.Y., Tamaki, S.J., and Keen, NT. (1989). Cloned avirulence genes from 

570 the tomato pathogen Pseudomonas syringae pv. tomato confer cultivar specificity 

571 on soybean. Proc Natl Acad Sci USA 86, 157-161 



Baltrus MPMI 



24 



Downloaded from http://biorxiv.org/on September 18, 2014 



572 Kvitko, B.H., Park, D.H., Velasquez, A.C., Wei, C.-F., Russell, A.B., Martin, G.B., 

573 Schneider, D.J., and Collmer, A. (2009). Deletions in the Repertoire of 

574 Pseudomonas syringae pv. tomato DC3000 Type III Secretion Effector Genes 

575 Reveal Functional Overlap among Effectors. PLoS Pathog 5, e1 000388. 

576 Lee, A.H., Hurley, B., Felsensteiner, C, Yea, C, Ckurshumova, W., Bartetzko, V., 

577 Wang, P.W., Quach, V., Lewis, J.D., Liu, Y.C., Bornke, F., Angers, S., Wilde, A., 

578 Guttman, D.S., and Desveaux, D. (2012) A bacterial acetyltransferase destroys 

579 plant microtubule networks and blocks secretion. PLoS Pathog 8, e1 002523 

580 Lemonnier, M., Landraud, L, and Lemichez, E. (2007). Rho GTPase-activating bacterial 

581 toxins: from bacterial virulence regulation to eukaryotic cell biology. FEMS 

582 Microbiol Rev 31, 515-534. 

583 Lindeberg, M., Cunnac, S., and Collmer, A. (2012). Pseudomonas syringae type III 

584 effector repertoires: last words in endless arguments. Trends Microbiol 20, 

585 199-208. 

586 Melotto, M., Underwood, W., Koczan, J., Nomura, K., and He, S.Y. (2006). Plant 

587 stomata function in innate immunity against bacterial invasion. Cell 126, 969-980. 

588 Morris, C.E., Sands, D.C., Vinatzer, B.A., Glaux, C, Guilbaud, C, Buffiere, A., Yan, S., 

589 Dominguez, H., and Thompson, B.M. (2008). The life history of the plant 

590 pathogen Pseudomonas syringae is linked to the water cycle. ISME J 2, 321 -334. 

591 Mucha, E., Fricke, I., Schaefer, A. Wittinghofer, A., and Berken, A. (201 1) Rho proteins 

592 of plants— functional cycle and regulation of cytoskeletal dynamics. Eur J Cell 

593 Biol 90, 934-943 

594 Mucyn TS, Yourstone S, Lind AL, Biswas S, Nishimura MT, Baltrus DA, Cumbie JS, 

595 Chang JH, Jones CD, Dangl JL, Grant SR (2014) "Variable suites of non-effector 

596 genes are co-regulated in the type III secretion virulence regulons across the 

597 Pseudomonas syringae phylogeny" PLoS Pathogens. 10, e1 003807. doi: 

598 10.1371/journal.ppat.1003807. 
599 

600 Nakamura, S., Mano, S., Tanaka, Y., Ohnishi, M., Nakamori, C, Araki, M., Niwa, T., 

601 Nishimura, M., Kaminaka, H., Nakagawa, T., Sato, Y., and Ishiguro, S. (2010). 

602 Gateway binary vectors with the bialaphos resistance gene, bar, as a selection 

603 marker for plant transformation. Biosci Biotechnol Biochem 74, 1315-1319 

604 

605 O'Brien, H.E., Thakur, S., and Guttman, D.S. (201 1). Evolution of Plant Pathogenesis in 

606 Pseudomonas syringae: A Genomics Perspective. Annu Rev Phytopathol 49, 

607 269-289. 

608 Quigley, N.B., and Gross, D.C. (1994). Syringomycin production among strains of 

609 Pseudomonas syringae pv. syringae: conservation of the syrB and syrD genes 



Baltrus MPMI 



25 



Downloaded from http://biorxiv.org/on September 18, 2014 



610 and activation of phytotoxin production by plant signal molecules. Mol Plant 

611 Microbe Interact 7, 78-90. 

612 Ramos, C, Matas, I.M., Bardaji, L, Aragon, I.M., and Murillo, J. (2012). Pseudomonas 

613 savastanoi pv. savastanoi: some like it knot. Mol Plant Pathol 13,998-1009 

614 Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D.L., Darling, A., Hohna, S., Larget, 

615 B., Liu, L, Suchard, M.A., and Huelsenbeck, J. P. (2012). MrBayes 3.2: Efficient 

616 Bayesian Phylogenetic Inference and Model Choice Across a Large Model 

617 Space. Systematic Biology 61, 539-542. 

618 Schellenberg, B., Ramel, C, and Dudler, R. (2010). Pseudomonas syringae virulence 

619 factor syringolin A counteracts stomatal immunity by proteasome inhibition. Mol 

620 Plant Microbe Interact 23, 1287-1293. 

621 Studholme, D.J. (201 1 ). Application of high-throughput genome sequencing to 

622 intrapathovar variation in Pseudomonas syringae. Mol Plant Pathol 12,829-838 

623 Vinatzer, B.A., Teitzel, G.M., Lee, M.-W., Jelenska, J., Hotton, S., Fairfax, K., Jenrette, 

624 J., and Greenberg, J.T. (2006). The type III effector repertoire of Pseudomonas 

625 syringae pv. syringae B728a and its role in survival and disease on host and 

626 non-host plants. Mol Microbiol 62, 26-44. 



627 

628 Acknowledgements 

629 Funding was provided by startup funds to David Baltrus from the University of Arizona. 

630 We thank the University of Delaware Sequencing Core Facility for technical help with 

631 PacBio sequences. We thank Cindy Morris, INRA, France for kindly providing access to 

632 strains described in this manuscript. We also thank Jeff Dangl for experimental reagents 

633 and positive encouragement. 



Baltrus MPMI 



26 



Table 1. Strains and Plasmids 



Strain 


Description 


Reference 


CC1466 






CC1557 


P. syringae CC1557 


Morris et al. 2008 


CC1583 


P. syringae CC1583 


Morris et al. 2008 


DBL627 


A rifampicin resistant version of P. syringae CC1557 


This paper 


DBL704 


DBL627 with pDAB334, CC1557 rifR mutant with EV 


This paper 


DBL690 


DBL627 with pDBL048 integrated, hopBJl deletion construct integrated into CC1557 rifR 


This paper 


DBL692 


DBL627 with pDBL051 integrated, hrcC deletion construct integrated into CC1557 rifR 




DBL698 


DBL627 hrcC mutant 


This paper 


DBL696 


DBL627 hopBJl mutant 


This paper 


DBL704 


DBL627 with EV pDAB334, CC1557 rifR with EV 


This paper 


DBL706 


DBL696 del with pDAB334, hopBJl deletion with EV 




DBL705 


DBL696 del with pDBL052, hopBJl deletion with hopBJl 


This paper 


DBL708 






DBL707 


DBL698 with pDBL052, CC1557 rifR hrcC- with hopBJl 


This paper 


DBL604 


PtoDC3000 with pDBL061 


This paper 


DBL888 


PtoDC3000 hrcC mutant with hopBJl (no promoter) pDBL061 


This paper 


MTN3101 


GV3101 with pMTN3100; 35S HopBJl:YFP 




MTN3131 


GV3101 with pMTN3100; 35S HopBJl:YFP 


This paper 


MTN3132 






MTN2890 


PtoDC3000 with pDBL052; HopBJl in pJC531 


This paper 


MTN3151 


PtoDC3000 with pMTN3149; HopBJl C174S in pJC531 


This paper 


MTN3152 


PtoDC3000 with pMTN3150; HopBJl H192A in pJC531 


This paper 


MTN3155 


DBL627 with pMTN3149; HopBJl C174S in pJC531 




MTN3156 


DBL627 with pMTN3150; HopBJl H192A in pJC531 


This paper 


MTN3240 






MTN3241 


PtoDC3000-D28E with pMTN3225; HopBJl in pJC532 


This paper 



r IdSmiU 


uescripiion 


MnilDIOLICS 


Keierence 


pMTN1907 


Gateway destination vector for making clean deletions in P. syringae 


TetR, KanR, CamR, 
SucS 


Baltrus et al., 2012 


pDONR207 


Gateway entry vector from Invitrogen 


GentR, CamR 


Invitrogen 


pMTN1970 


the 3' UTR of Arabidopsis HSP18.2 in pDONR207 


GentR 


Baltrus etal., 2012 


pJC531 


Gateway destination vector based off of pBBRMCSl no promoter; C-terminal HA tag 


KanR, CamR 


Chang et al., 2005 


pJC532 


Gateway destination vector based off of pBBRMCSl no promoter; C-terminal AAvrRpt2 


KanR, CamR 


Chang et al., 2005 


pDAB334 


MTN1970-2 EV construct in pJC531 


KanR 


Baltrus et al., 2012 


pDBL040 


hopBJl cloned with promoter, no stop codon in pDONR207 


GentR 


This paper 


pMTN3111 


hopBJl C174S with promoter, no stop codon in pDONR207 


GentR 


This paper 


pMTN3112 


hopBJl H192A with promoter, no stop codon in pDONR207 


GentR 


This paper 


pDBL052 


hopBJl cloned with promoter, no stop codon in pJC531 


KanR 


This paper 


pMTN3149 






This paper 


pMTN3150 


hopBJl H192A cloned with promoter, no stop codon in pJC531 


KanR 


This paper 


pDBL042 


hopBJl deletion construct in pDONR207 


GentR 


This paper 


pDBL048 


hopBJl deletion construct in MTN1907 


TetR 


This paper 


pDBL012 


hopBJl cloned without promoter into pDONR207, no stop 


GentR 


This paper 


pMTN3113 


hopBJl C174S ORF, no stop codon in pDONR207 


GentR 


This paper 


pMTN3114 






This paper 


pDBL061 


hopBJl cloned without promoter into pBAV187 


TetR 


This paper 


pMTN3224 


HSP18.2 UTR in pJC532 


KanR 


This paper 


pivl 1 INoZZj 


nopDJi cionea wiin promoter, no stop coaon in pj^ooz 




i nis paper 


pMTN3226 


hopBJl C174S cloned with promoter, no stop codon in pJC532 


KanR 


This paper 


pMTN3227 


hopBJl H192A cloned with promoter, no stop codon in pJC532 


KanR 


This paper 


pDBL045 






This paper 


pDBL051 


CC1557 hrcC deletion construct in MTN1907 


TetR 


This paper 


pGWB641 






Nakamura et al., 
2010 


pMTN3100 


hopBJl ORF in pGWB641 


SpecR 


This paper 


pMTN3119 


hopBJl C174S ORF in pGWB641 


SpecR 


This paper 
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635 



Table 2. Type III Effector Family Distribution 



avrE 
hopF2-l 
hopF2-2 
hopMl 
hopYl 
hopAAl 
hopAGl 
hopAHl 
hopAH2 
hopAll 
hopASl 
hopBFl 
hopBJl 



CC1466 




636 



The letter P denotes presence of this particular effector family, whereas the lack of a letter 
denotes absence. 



637 
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Table 3. Responses of Arabidopsis Accessions to HopBJl 



Arabidopsis Accession 



Reaction 




Kn-0 

Ler-0 

Mt-0 

No-0 

Oy-0 

Po-0 

Rsch 

Sf-0 

Tsu-0 

Wil-2 

Ws-0 

Wu-0 

Zu-0 




639 



++ HR in 24hr 
+ HR in 48hr 

Gene names next to accession 
indicate mutant lines 



640 
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641 

642 Figure 1. P. syringae CC1557 Virulence in N. benthamiana is Type III Secretion 

643 Dependent. A) Growth of wild type P. syringae CC1 557 as well as a hrcC mutant was 

644 measured in planta three days post syringe inoculation. Measurements are based on 

645 three independent experiments with at least 4 replicates in each. Bacterial population 

646 sizes are significantly different (Tukey's HSD, p<0.05) between strains. Error bars 

647 represent 2 standard errors. Dashed lines indicate approximate population sizes at day 

648 0, which were not significantly different. 
649 

650 Figure 2. Diversity of the CEL locus across P. syringae sensu latu. A) Bayesian 

651 phylogenies for AvrE (left), HrpW (middle), and concatenated MLSA loci (right), for a 

652 subset of P. syringae strains. Posterior probabilities for all nodes are >0.95. Scale bars 

653 indicate number of amino acid changes. Phylogenetic patterns for both of these loci 

654 approximate relationships based on core genomes and MLSA loci, except that HrpW 

655 from P. syringae CC1583 clusters more closely with P. syringae pv. tomato DC3000. 

656 This suggests a recombination event at this locus within P. syringae CC1583. B) 

657 Genomic context of the CEL across P. syringae sensu latu. In most strains, hopM1 is 

658 bordered by avrE and hrpW. Some P. viridiflava strains lack hopM1 (Bartoli et al., 2014), 

659 while CC1557 and CC1466 have replaced this region with hopBJL Arrows in represent 

660 ORFs while the gap and triangle represent a deletion of this region within strains 

661 CC1524 and CC1417. 
662 
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663 Figure 3. HopBJI is a Virulence Factor for P. syringae CC1557 in N. benthamiana. 

664 We created a deletion mutant of hopBJI in P. syringae CC1557, and complemented this 

665 strain by expressing hopBJI from a plasmid using the native promoter. Deletion of 

666 hopBJI lowers growth of P. syringae CC1557 three days post inoculation. Both of these 

667 phenotypes are complemented by expression of wild type hopBJI in trans. Three 

668 independent experiments were performed with at least 4 replicates in each treatment. 

669 Letters represent statistical differences within an ANOVA (Tukey's HSD, p<0.05). Error 

670 bars represent 2 standard errors. Dashed lines indicate approximate population sizes at 

671 day 0, which were not significantly different. B) Leaves from two-week old N. 

672 Benthamiana plants were syringe inoculated with various strains derived from P. 

673 syringae CC1 557. Three days after inoculation, plants showed evidence of tissue 

674 collapse dependent on a functioning T3SS (compare wild type vs. hrcC mutant) as well 

675 as HopBJI (compare hopBJI deletion strain vs. hopBJI complement). Picture is 

676 representative of multiple trials. 
677 

678 Figure 4. Protein Motifs and catalytic residues within HopBJI are similar to the E. 

679 coll CNF1 Toxin. A) Despite no sequence similarity between HopBJI and CNF1 by 

680 BLASTp searches, modeling of tertiary protein structure demonstrates limited similarity 

681 in protein structure between the two. Importantly, the catalytic dyad of E. coli CNF1 

682 (Cys866-His881) perfectly lines up with a putative catalytic dyad within HopBJI (Cys174 

683 and His192). B) All constructs of HopBJI can be translocated by Pto DC3000 D28E 

684 (Cunnac et al., 201 1) into Arabidopsis. Top: translocation tests into accession Col-0. 
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685 Bottom: translocation tests into rps2-101C mutant plants. Translocation of AvrB is 

686 shown as a positive control and causes HR in both plant backgrounds, while the Pto 

687 DC3000 D28E mutant causes no cell death. HopBJI causes rapid cell death across 

688 both plant backrounds when translationally fused to the C-terminus of AvrRpt2, while 

689 catalytic mutants only show cell death in Col-0. Picture is representative of >4 

690 independent leaves per strain / plant combination C) Catalytic mutants of HopBJI 

691 cannot complement the growth defect of a CC1 557 hopBJI deletion strain. D) 

692 Translocation of HopBJI , but not either of the two catalytic mutants, lowers growth of 

693 Pto DC3000 on Arabidopsis Col-0. E) Wild type and the complemented AhopBJI 

694 deletion of strain CC1 557 cause tissue collapse in N. benthamiana. Neither of the 

695 catalytic mutants complements the loss of virulence of the AhopBJI deletion strain. 

696 Picture is representative of multiple assays. F) Transient expression of HopBJI , but not 

697 the catalytic mutants, by Agrobacterium in N. benthamiana causes cell death. Western 

698 blots and Ponceau staining for each of these Agrobacterium infiltrations are displayed 

699 below. Picture is representative of multiple assays. Each growth curve experiment 

700 consisted of at least two independent trials, but representative data from only one trial is 

701 shown in each case. Letters indicate significant differences between bacterial population 

702 sizes between strains (p<0.05, Tukey's HSD). Dashed lines indicate approximate 

703 population sizes at day 0, which were not significantly different. 
704 

705 Figure S1. Negative Correlation Between Type III Effector Number and Presence 

706 of Syringolin Pathways. The number of type III effectors is plotted for a subset of 



Baltrus MPMI 



35 



Downloaded from http://biorxiv.org/on September 18, 2014 



707 strains investigated within Baltrus et al. 201 1 as well as CC1557 and CC1583, with 

708 strains grouped by the presence of a syringomycin-like pathway. Full strain names are 

709 described in Baltrus et al. 201 1 . Genomes with the genetic potential to encode a 

710 syringomycin-like toxin contain a significantly lower number of type III effectors 

711 (Wilcoxian rank sum test, p<0.001 5). 
712 

713 Figure S2. P. syringae CC1557 grows as well as P. syringae pv. syringae B728a 

714 on N. benthamiana. Bacterial population sizes after one day of growth are shown for 

715 wild type and hrcC- versions of the established N. benthamiana pathogen P. syringae 

716 pv. syringae B728a. Also shown are population sizes for wild type and hrcC- versions of 

717 strain CC1 557 as well as wild type P. syringae pv. tomato DC3000. Dashed line 

718 indicates approximate day 0 population sizes for all strains, which were not significantly 

719 different than one another. Error bars represent 2 standard errors. Growth curves were 

720 performed twice with similar results, but only data from one trial are shown. Letters 

721 indicate significant differences in population sizes between strains (p<0.05, Tukey's 

722 HSD) 
723 

724 Figure S3. Toxin Regions are Syntenic in P. syringae pv. syringae B728a and P. 

725 syringae CC1557. Mauve was used to align four whole P. syringae genomes, and the 

726 region that contains pathways encoding syringomycin was extracted from P. syringae 

727 pv. syringae B728a. The toxin region for syringomycin and syringopeptin is shown in 

728 purple, while adjacent genomic regions conserved throughout the four genomes are 
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729 listed with either the number 1 or 2. The strain order is (from top): CC1557, pv. syringae 

730 B728a, pv. phaseolicola 1448a, pv. tomato DC3000 

731 

732 Figure S4. Putative hrp-box for hopBJI in P. syringae CC1557. The sequence of the 

733 coding strand around hopBJI is shown from strain CC1 557. The ORF of hopBJI is 

734 shown as an arrow. The region immediately upstream the hopBJI start codon is also 

735 shown, with the predicted hrp-box and consensus nucleotides labeled in red. 

736 
737 
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